The delayed uncoupled continuous-time random walks 
do not provide a model for the telegraph equation 
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It has been aUeged in several papers that the so caUed delayed continuous-time random walks 
(DCTRWs) provide a model for the one-dimensional telegraph equation at microscopic level. This 
conclusion, being widespread now, is strange, since the telegraph equation describes phenomena with 
finite propagation speed, while the velocity of the motion of particles in the DCTRWs is infinite. 
In this paper we investigate how accurate are the approximations to the DCTRWs provided by the 
telegraph equation. We show that the diffusion equation, being the correct limit of the DCTRWs, 
gives better approximations in L2 norm to the DCTRWs than the telegraph equation. We conclude 
therefore that, first, the DCTRWs do not provide any correct microscopic interpretation of the one- 
dimensional telegraph equation, and second, the kinetic (exact) model of the telegraph equation is 
different from the model based on the DCTRWs. 

PACS numbers: 05.40.-a, 05.40.Fb, 05.60.-k, 05.60.Cd, 02.30.Jr, 02.30.Mv 



I. INTRODUCTION 

The continuous-time random walks (CTRWs) [ll-0| 
present sufficiently wide and general class of random 
walks, describing random motions of particles, when both 
the waiting (sojourn) time between successive jumps and 
the jump length (including its direction) are generally 
coupled random variables. The CTRWs have been widely 
applied in building up models of anomalous diffusion and 
transport in physics [sl-fisj (to name a few) and biology 
0, [3 [IE] , and in economical problems ^Td\ . 

In their original form ^j] the CTRWs describe random 
walks, or, more precisely, jumps, on lattices when only 
the waiting time is a random variable. In the general case 
of the uncoupled CTRWs both waiting time and jump 
length are independent random variables, this leads to 
the jump model of the CTRWs. The kinetic model of ran- 
dom walks, often called random flights, originates from 
the Pearson-Rayleigh random walks [4j. In this model 
the particle moves with constant velocity for a random 
time along straight lines between points where it changes 
randomly the direction of movement. The random flights 
can be considered both in the framework of the general- 
ized linear transport (kinetic) equation jl7| and in the 
framework of the CTRWs, when the waiting time and 
jump length are coupled random variables. The latter is 
the velocity or velocity-jump model of the CTRWs [Tsj . 

The jump model of the CTRWs has gained high pop- 
ularity due to its greater simplicity than the kinetic 
(velocity-jump) model. However, in the jump model the 
waiting time can generally be arbitrarily small or/and 
the jump length can be arbitrary long (Levy flights), and, 
therefore, a walking particle can move with infinite veloc- 
ity. Thus, the jump model may, in contrast to the kinetic 
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one, violate the principle of causality. 

This is similar to the classic diffusion paradox [l9l - [23j 
related to the conventional diffusion equation, which is 
widely used for approximate macroscopic description of 
nonanomalous diffusion and Brownian motion. The dif- 
ference is that in classic Einstein's model of Brownian 
motion ^24|, leading to the diffusion equation, the parti- 
cle moves with finite velocity. In fact Einstein's model 
is nothing but the symmetric Bernoulli random walk 
[2^ , which can be considered as a degenerate case of the 
CTRWs in the frameworks of both the jump and kinetic 
models. 

The telegraph equation [H [13, [H [H Ssl was pro- 
posed as an alternative to the diffusion equation. In 
contrast to the latter, which is parabolic, the telegraph 
equation is hyperbolic, providing the finite speed of sig- 
nal propagation. Two- and three-dimensional telegraph 
equations meet some formal problems since Green's func- 
tions can become negative [23|- Though in the case of 
thermal conduction this can be fixed by imposing re- 
strictions on the heat flux [2^. Yet the one-dimensional 
telegraph equation avoids the diffusion paradox and pro- 
vides better model for nonanomalous diffusion than the 
one-dimensional diffusion equation. 

In Ref. [2^ the authors discussed microscopic mod- 
els of the telegraph equation. Aside from the kinetic 
model they proposed another model based on uncoupled 
CTRWs with the waiting time distributed according to 
the gamma law {t/O"^) e~*/^ and the jump length hav- 
ing the finite second moment. The authors called the 
CTRWs with the gamma-distributed waiting time de- 
layed CTRWs (DCTRWs). Similar arguments were used 
in Ref. [1^ when deriving the telegraph equation with re- 
action. In Ref. [SO^I the telegraph equation with reaction 
was derived formally from the CTRW with exponentially 
distributed waiting time, and later, in Ref. [3ll |. — from 
the CTRW with a general waiting time. 

Earlier, the telegraph equation was derived in Ref. (32j 
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from uncoupled CTRWs with the waiting time having 
more general distribution than the gamma law. We 
call these CTRWs for brevity also delayed CTRWs (DC- 
TRWs). 

The conclusions that the DCTRWs simulates the tele- 
graph equation at the microscopic level are strange. In- 
deed, the telegraph equation describes phenomena with 
finite propagation speed, while the particle in the DC- 
TRWs may move with infinite velocity, i. e., the DCTRWs 
do not refiect the principal peculiarity of the telegraph 
equation. In spite of this discrepancy, one may believe 
that the telegraph equation still gives a precise approxi- 
mation to the DCTRWs. 

In this paper we investigate how accurate are the ap- 
proximations to the DCTRWs provided by the telegraph 
equation. We show that the diffusion equation, being the 
correct limit of the DCTRWs, gives better approxima- 
tions to the DCTRWs than the telegraph equation. One 
can conclude, therefore, that the DCTRWs cannot be a 
model for the telegraph equation at microscopic level. 

The rest of the paper is organized as follows. In 
Section |IT] we recall some facts concerning the uncou- 
pled CTRWs. Section [ITT] shows how the telegraph 
equation was derived from the DCTRWs in Ref. [2^ 
|32 |. The asymptotic behaviour of the DCTRWs is dis- 
cussed in Section |lVl In Section |V] we compare the tele- 
graph and diffusion approximations to two (discrete- 
and continuous-space) DCTRWs, completely described 
in Section \^ numerically, while asymptotic comparison 
of the approximations to the continuous-space DCTRW 
at long times is performed in Section |B] Section |VI] con- 
tains some conclusive remarks. 



II. UNCOUPLED CONTINUOUS-TIME 
RANDOM WALKS 



In this section we briefly formulate some facts concern- 
ing the CTRWs for further reference. 

We adopt the one-dimensional model and consider un- 
coupled CTRWs, when the waiting time and jump length 
are independent random variables, i.e., the joint proba- 
bility density if{x, t) for a particle to jump a distance x 
after waiting a time t is given by 



density pp^(x, t) is given by the integral equation |10l.l35 



ip{x,t) = ip{t)X{x), X e 



t > 0, 



where ■0(t) and A(x) are probability densities for the wait- 
ing time and jump length, respectively, with obvious nor- 
malizations /g 'ip{t) dt — 1 and ^i^) = 1. 

We consider here CTRWs in a continuous space, how- 
ever, they can be similarly described on a lattice, as it was 
originally proposed in pj], with obvious modifications. 

The CTRW is described by the density 
which is a probability density that the particle is at the 
point X at the time t, so that p^yj^{x, t) da; = 1. The 



i^t poo 



ip{x — x' ,t — t')p^^jif{x\t') dx' dt' 



W{t)5{x) 



+ / ^{t-t') 
Jo 



X{x - x')py^{x' , t') dx' 



dt', 
(II.l) 



where S{-) is the Dirac delta function, and 



tP{t')dt 



ij{t')dt, t>0, 



is the survival probability, i.e., the probability that a 
particle stays at the same position for the time t. Here 
one assumes that at the initial moment < = the particle 
was at the point a; = 0, i. e., Pp[^|t=o = S{x). 

The Fourier-Laplace transform of the equation pi.lt . 
denoted by J^C-, implies the MontroU- Weiss formula [l|. 



^Lp^^{k, s)- ^_ ^^^^^^ s)~ s[l~ FCv{k, s) 



1 - Cij{s) 



s[l-T\{k)Li){s)] 

where 

/oo 
\{x) e"=^ da; 
-OO 

is the Fourier transform of A, 

poo 

£^{s) = / i/.(i)e-"*di 



(II.2) 



is the Laplace transform of ip. 

In the uncoupled CTRWs the density p^^ can be ex- 
pressed explicitly (such a CTRW is the random walk sub- 
ordinated to a renewal process) [H, HI] : 



n=0 

oo 

= ^{t)S{x) +J2Pn{t)>^n{x) (II.3) 
n=l 

(ppj^(-, t), with fixed t, is actually a probability density of 
a random sum), where Pn(i) is the probability of n jumps 
occurring up to the time i, and A„(x) is the probability 
density of a distance x from the initial position reached by 
a particle after n jumps. The probability density An (a;) 
is given by 

A„(x) = A"*(a;), 
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where * means the convolution, A"* means the n-fold 
convolution of A with itself, i.e., X'^*{x) = S{x), X^*{x) = 
\{x), and A"*(a;) = (A("-i)* * A)(x). Equivalently, 

^A„(fc) - [J-A(fc)]". 

The probability Pn (t) is given by 

(clearly ip{t) — and '^{t) — for t < 0), in particular, 
PQ(t) ^^{t). Equivalently, 

£p^{s)^C^{s)[C^j{s)r. 

This completes the set of auxiliary propositions neces- 
sary for further considerations. 



III. TELEGRAPH APPROXIMATIONS TO THE 
DELAYED CONTINUOUS-TIME RANDOM 
WALKS 

Consider a family of probability densities for waiting 
time t m, m, [Hi 



e-*/^sinh ( , 0<a<l, 

a = 1, 

(III.l) 



(obviously, i/'il^) = linia-j-i V'al^))- The Laplace trans- 
form of the density i\}a is 



a 



93 + ly - {1 - a) 
a 



612 (s + l/6')2- (l-a)/6'2' 

The density was used in Refs. [H, (see also 
Ref . (S^l ) , and it belongs to the family of the gamma dis- 
tributions. The density -01 is a probability density of a 
sum of two independent exponentially distributed ran- 
dom variables with the probability density (1/0) e~*/^, 
and CTRWs with such a distribution of waiting times 
were called in Ref. ^ the delayed CTRWs (DCTRWs). 
Note that the mean waiting time for the exponential dis- 
tribution is 9, and it is 29 for the distribution -01. The 
mean waiting time for the distribution tpa is 20/a, i.e., 
the mean waiting time tends to infinity as a — >■ 0. Com- 
pletely understanding the conditional character of this 
notion we call for brevity the CTRWs with the probabil- 
ity densities V'a of the waiting time also the DCTRWs. 
We denote densities describing the DCTRWs by Pdr^. 
The Montroll- Weiss formula (jn.2[) with Cipa implies 
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The straightforward calculations show that the density 
Pdrw ^ solution to the equation 

2 dt^ dt 

OO 

A(x - x')p^^^{x' ,t) dx' - PDRw(a^'0 

-OO 

(III.2) 



a 
29 



with the initial value conditions 



PDRwlt=0 ~ '5 (2^) I 



dp 



DRW 



dt 



= 0. 



Consider a random walk on a regular one-dimensional 
lattice, where the distance between nearest-neighbour 
points is (T, and the coordinates of the lattice points are 
xi = la, I ^ Z. Suppose that jumps are made to the 
nearest-neighbour points with equal probability. This 
discrete-space random walk is described by the proba- 
bility density Ad (a:) = [5{x — a) + 5{x + (j)]/2, and the 
equation pn.2[) becomes 

^^Vdrw(^^^) , dp^^^{xi,t) 



dt 



a 
29 



PDRwi^l - CT, + PuRwi^l + 



(in.3) 



Assuming that the density p-^-^w smooth with respect 
to x, one can expand it in the Taylor series in a;. If 
a is small enough then we neglect the derivatives of an 
order higher than two, and obtain the approximation for 
the density Pqj^^ given by the solution of the telegraph 
equation (TE) 



9_ S^Pte 
2" 



dp. 



TE 



dt^ dt 
with the initial conditions 



~A9 dx^ 







PTE I 



Six), 



dp. 



TE 



dt 



0. 



(IIL4) 



(IIL5) 



t=o 



s{9s 



a[F\{k) - I]/ 9' 



The same result is obtained if jumps are made not only to 
the nearest-neighbour points, the only decisive condition 
is that the variance of the distribution Ad is to be equal 
to 0-2. 

Similar derivation of the telegraph equation from the 
discrete-space random walk with the probability density 
il>a for waiting time was performed in Ref. [32j , where the 
authors wrote about "the continuum limit" as cr — > 0. 

Note that both the equations pL2|) and pL3|) de- 
scribe particles, moving with infinite velocity, while the 
telegraph equation describes phenomena with finite prop- 
agation speed. 

Another derivation of the telegraph approximation is 
given in Ref. psj . This is similar to the above one, but it 



4 



is performed in the Fourier-Laplace space. The deriva- 
tion, given in Ref. [l^ for a = 1, is in fact as follows. If 
the distribution of the jump length is symmetric and has 
the finite second moment with the variance equal to a^, 
then the Fourier transform of the probability density of 
the jump length is approximately J-X{k) « 1 — for 
\ak\ <C 1. After substitution of this J^A(fc) and C4>a{s) 
into the MontroU-Weiss formula (jll.2[) one obtains that 
the Fourier-Laplace transform of the density Pqj^^ ^P" 
proximately equal to 



J^^PuRwi^^ ~ ^^PT^{k, s) for \ak\ < 1 



with 



(in. I 



The inverse Fourier-Laplace transform of TLprj--^ implies 
that the density p^^^ is a solution to the telegraph equa- 
tion (jIII.4[) with the initial conditions (IIII.5I) . The au- 
thors of the paper call the telegraph equation pil.4[) 
"the limit of small jumps (diffusive approximation)" of 
the CTRW with the probability density ipi (jlll.ip for 
waiting time. 

The conclusion that the solution to the problem (IIII.4I) , 
piI.Sp approximates the DCTRW is strange, since the 
telegraph equation describes phenomena with finite prop- 
agation speed. At the same time, the velocity of the 
motion of the particle in the DCTRW is infinite, since 
waiting time can be arbitrarily small. Moreover, the 
telegraph equation obviously is not "the limit of small 
jumps" of the DCTRWs. 



IV. THE ASYMPTOTIC BEHAVIOUR OF THE 
DELAYED CONTINUOUS-TIME RANDOM 
WALKS 

The asymptotic behaviour of various CTRWs was 
studied in Refs. [H, [s^ H^- For convenience we briefly 
repeat the as ymp totic analysis for the DCTRWs on the 
basis of Refs. [s^ S|- To find the asymptotic behaviour 
of the DCTRWs the densities X{x) and V'a(^) ^'^^ replaced 
by the scaled densities \h{x) — \{x/h)/h, h > 0, and 
i'a,T{t) = il)g_{t/T)/T, T > 0, respectively [H, [s^. The 
parameters h and r can be considered as the charac- 
teristic step length and waiting time, respectively. The 
Fourier transform of \h is F\h{k) — J-X{hk), and the 
Laplace transform of -00. r is £V'a,r(s) = Cil^^iTs). We 
suppose that the distribution of the jump length is sym- 
metric, A(— a;) = A(x), and has the finite second moment 
with the variance equal to cr^, then the asymptotic be- 
haviour for J^Xh is 1 — TXh{k) ^ h'^a'^k'^/2 as /i — ^ 0, 
fc e M [39]. For the density tpa,T one has the asymptotics 
1 — Cipa^ris) ^ 2T9s/a as r — 7^ 0. The latter asymp- 
totics are also valid in general case for V'r(0 = '4>{t/T)/T, 
if the density V' has the finite expectation equal to 26'/a 
[39I . After substitution of the above asymptotics into the 



Montroll- Weiss formula (jll.2p . one can see that the only 
possible nontrivial limit is 



1 - C'lpa,Tis) 



s[l~TXh{k)C^aAs)] 
— !• TCpY)^{k, s) as /i and r ^ (IV.l) 

under the scaling relation h'^/r — 1, where 

1 



J"£pDE(fc, s) 



72fc2/(46') 



The inverse Fourier-Laplace transform implies that the 
density p^g is a solution of the diffusion equation (DE) 



dt AO dx^ 
with the initial condition 



= 



PdeI 



6{x) 



(IV.2) 



(IV.3) 



Thus, the asymptotic behaviour of the DCTRWs is 
given by the diffusion equation (jIV.2p . Note that the dif- 
fusion equation is obtained from the telegraph equation 
(jIII.4p by omitting the second time derivative. 

Note also that J^£prj.^{0, s) — 1/s and J"£pqe(0, s) — 
1/s, or, equivalently, J^prpg(0,t) — 1 and TpY^^{0,t) = 1 
for t > 0, i. e., the law of conservation of particles is valid 
for both the telegraph and diffusion approximations. 

The other possible limits are TCpj^^ f^^{k,s) — > 1/s 
as ft. — ?> 0, and TCpj^^ ^ r(^j'^) ~^ 1/* as r — 7> and 
h/y^ — > 0. The former is the limit of infinitesimal jumps 
with finite mean waiting time. This is the correct limit 
[to be] obtained in Refs. [H,!!^- These asymptotics are 
trivial, since the inverse Fourier-Laplace transform of 1 /s 
gives S(x), t > 0, which means that the particle does not 
leave an initial position. 

The diffusion limit, under the scaling relation h'^ /t = 
const, is valid for a variety of symmetric (unbiased) ran- 
dom walks, both continuous-time and discrete-time ones 
0; [H; [13: • Ih partlcular, it is valid for the symmetric 
Bernoulli random walk, which is discrete both in time 
and space, when the waiting time is exactly r and the 
jumps are exactly ±h with equal probabilities, i.e., the 
corresponding probability densities are 

tpdit) = S{t - t), t > 0, and 

Xd{x) = ^[S{x - h) + S{x + h)], h>0. (IV.4) 

However, in contrast to the CTRWs, in which waiting 
time can be arbitrarily small and, hence, the particle 
moves with an infinite speed, the speed of the motion 
of the particle in the Bernoulli random walk is finite. 
Besides, it turns out that the solution of the telegraph 
equation approximates the density for the Bernoulli ran- 
dom walk better than that of the diffusion equation ■ 
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It is necessary to note here that the explanation, given in 
Ref. [iOl, does not allow us to judge whether the telegraph 
equation approximates the random walk better than the 
diffusion equation in a wide time interval? To check this 
conclusion we have performed calculations, in which the 
binomial distribution, being the distribution of the walk- 
ing particle, was evaluated through the first two terms of 
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Stirling's asymptotic series for the gamma function 
r(z) = v^27r/z(z/e)^ [l + l/(12z) + 0(|z|-2)] as \z\ oo, 
z G C, \ argz| < TT — A, A > 0, which is more accurate 
than conventional Stirling's formula. These calculations 
have confirmed the conclusion up to time t « lOOr, which 
seems to be asymptotic. Obviously, the solutions of both 
the diffusion equation and the telegraph one do not ap- 
proximate the binomial distribution at short time. 

The reason, by which the solution of the telegraph 
equation approximates the Bernoulli random walk bet- 
ter than that of the diffusion equation, seems to be clear: 
the Bernoulli particle moves with finite velocity. How- 
ever, the probability density for the waiting time ip^, 
Eq. (jIV.4[) , can be weakly approximated with arbitrary 
accuracy by the gamma distribution: 



1 



fit 

T 

weakly 



> S{t - t) as ^ +00, (IV.5) 

where T is the Gamma function; this can easily be derived 
with the Laplace transform. In the random walk with 
the distribution -0^ of the waiting time the particle moves 
with infinite velocity, since the waiting time can generally 
be arbitrarily small. However, for sufficiently large fi the 
telegraph approximation to this random walk is clearly 
better than the diffusion one, since ipfj.it) is very "close" 

to S(t~T). 

Thus, in spite of the diffusion asymptotic behaviour 
for the DCTRWs and the infinite speed of the motion 
of the particle, the question still remains: which of the 
telegraph or the diffusion equations gives better approx- 
imation to the DCTRWs? 



V. COMPARISON OF THE TELEGRAPH AND 
DIFFUSION APPROXIMATIONS TO THE 
DCTRWS 

To compare the telegraph and diffusion equations with 
the DCTRWs we consider two particular DCTRWs. One 
of them is a discrete-space random walk, the other is the 
continuous-space one. 

The discrete- space random walk takes place on a reg- 
ular one-dimensional lattice, where the distance between 
nearest-neighbour points is a, and the coordinates of the 
lattice points are xi = la, I £ Z. Jumps are made to 
the nearest-neighbour points with equal probability. This 
random walk is described by the probability density 

Mx)^\[6{x-a) + 6{x + a)] . (V. 1) 



The subscript b means Bernoulli (since this is the same 
distribution as in the Bernoulli random walk). The vari- 
ance of this distribution is equal to a^. 

The continuous- space random walk has the Gaussian 
distribution of the jump length with the variance equal 
also to a^: 



AG(a;) 



1 



/27r(T 



(V.2) 



its Fourier transform is J-XQ{k) 



''^'Z^. The Gaus- 



sian distribution of the jump length was used in Ref. [29| . 

The densities of the discrete- and continuous-space 
DCTRWs are given in Appendix E] by Eqs. (|A.1|) and 
(|A.2p . respectively. The telegraph and diffusion approx- 
imations are given by Eqs. (IA.5I) and (IA.8I) . respectively. 

All numerical results are obtained with the parameters 
a — 1 and 9 = 1. 

Figs.[T] and [5] show (in Cartesian and logarithmic 
scales) the densities of the continuous- and discrete-space 
DCTRWs and their telegraph and diffusion approxima- 
tions at intermediate values of time t = 5 and t = 10 with 
a = 1. Figs.[3]and|3]show differences between each of the 
two approximations and the densities of the continuous- 
and discrete-space DCTRWs, respectively, obtained for 
the same values of time and the parameters. Calculations 
show that in all the cases the L2-discrepancy for the dif- 
fusion approximation is less than that for the telegraph 
one. In the case of the discrete-space random walk we cal- 
culated the ^2-discrepancy. Note also that the maximum 
absolute value of the difference for the diffusion approxi- 
mation is in all the cases less than that for the telegraph 
approximation. Besides, it is important to emphasize 
that the telegraph approximation is incorrect with re- 
spect to the velocity of the motion of the particle, see 
Figs.[T] and [2] for i = 5. The telegraph model gives the 
finite velocity, while it is infinite in the DCTRWs. The 
diffusion model is correct in this respect, however the 
diffusion paradox remains. 




FIG. 1. (Color online) The densities Pnnwci^^'t) ^-nd 
PoRWBi^ht) and the telegraph and diffusion approximations 
PTE(2;,i) and pDE(a;,i) at t = 5 and t = 10 with a = 1, a = 1 
and 9=1. The singular terms Pdrwg ^-^d PxE ^re depicted 
by vertical arrows. 
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FIG. 2. (Color online) The same as in Fig.[Tl with logarithmic 
scale on the vertical axis. 
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FIG. 3. (Color onhne) The differences p\<^(x, t) — Pdrwg(^i ^) 
and pDE(a;,t) — Pdrwg(3^i*) at t = 5 and t = 10 with a — 1, 
a — 1 and 6 = 1. Note that the vertical scale unit is 10"^. 
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Fig. [5] shows graphs resulting from the Fourier trans- 
forms -T^PdrwgC^' -^PteC^' ^i^'i -^Pde(^' ^ = 20 
and t = 50 with a = 1. These values correspond to the 
values of the small parameter e — 0.05 and e = 0.02, re- 
spectively (see Appendix [BJ . Fig. |6] shows the integrands 
of the integrals (|B.ip for the same values of time and 
parameters a, a and 6 as in Fig.[SJ Note that the asymp- 
totics (IB.6|) and (IB.7|) provide very good approximations 
(not shown in the figure) to the integrands. 
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FIG. 5. (Color online) Graphs resulting from the Fourier 
transforms .7^Pdrwg(^i -^PteC^j^) ^nd J-'pj^^{k,t) a.t t = 
20 and t = 50 with a = 1, a = 1 and ^ = 1, i. e., e = 0.05 and 
£ = 0.02. 
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FIG. 6. (Color online) The integrands of the integrals (|B.1[) 
at t = 20 and t = 50 with cr = 1, a = 1 and 9 = 1, i.e., 
£ = 0.05 and e = 0.02. Note that the vertical scale unit is 
10-*. 



FIG. 4. (Color online) The differences pxE(2;!,i) 
Pdrwb (a;i,i) and pDE(a;;, t) -Pdrwb (a^!, i) at f = 5 and f = 10 
with a = 1, a = 1 and d = 1. Note that the vertical scale unit 
is 10"^ 



To assess the approximations at long times we 
have performed in Appendix |B] asymptotic comparison 
of the telegraph and diffusion approximations to the 
continuous-space DCTRW through their Fourier trans- 
forms. 



For comparison we show in Fig. [7] differences between 
each of the two approximations and the densities of 
the discrete-space DCTRW (not Fourier transforms), ob- 
tained at the same long times t — 2Q and t — with 
a = 1. Note that Fig. [7] is qualitatively similar to Fig. 21 
In these cases the Z2-discrepancy for the diffusion approx- 
imation is also less than that for the telegraph one. We 
do not show differences between the two approximations 
and the densities of the continuous-space DCTRW for 
i = 20 and t = 50 because Fig. IH] illustrates this via the 
Fourier transforms. 
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FIG. 7. (Color online) The differences p^^{xi,t) — 
PDBWBi^h't) and pj^js{xi,t) - PoBWEi^ht) at t ^ 20 and 
t = 50, a = 1, a = 1 and 6 — 1. Note that the vertical 
scale unit is 10"^. Cf. with Fig. ID 



Figs. [8] and [9] show the densities of the continuous- 
and discrete-space DCTRWs and their telegraph and 
diffusion approximations with the parameter a = 0.1. 
Fig.|S] corresponds to the small value of time t = 20 
(e — 0.5), which is equal to the mean waiting time. In 
this case the singular term Pdrvvg is quite "heavy" : the 
probability that the particle stays at the initial point is 
!f'o.i(20) ~ 0.37. Note that the telegraph and diffusion 
"approximations" essentially differ from the densities of 
the DCTRWs. At the same time they are very close 
to each other. Fig. [9] corresponds to intermediate and 
long time: t = 50, t ~ 100 (intermediate) and t — 500 
(long) (e = 0.2, e = 0.1 and e = 0.02, respectively). For 
t = 50 the singular term Pdrwg ^^^^^ "heavy" enough: 
the probability that the particle stays at the initial point 
for the time t = 50 is !P'o.i(50) w 0.08. It is necessary 
to emphasize that all the results are in agreement with 
the asymptotics (jB.8|) - (jB.10P (even for intermediate and 
small values of time) : the differences between the two ap- 
proximations are much less than the differences between 
each of them, and the densities of the DCTRWs, while 
the diffusion approximation is slightly better than the 
telegraph one. 

The DCTRWs with small a can be explained as fol- 
lows. Small values of a mean that the particle mostly 
remains in rest than jumps, since the mean waiting time 
is long enough. Small values of time can be described 
as lesser than or comparable with the mean waiting time 
26/a, i.e., such that at/d < 1. Therefore, for small val- 
ues of time both the diffusion and telegraph equations 
badly approximate the DCTRW, see Fig. [51 At the same 
time, if at/9 ~ 1 then t/9 ^ 1, i.e., the singular term 
(|A.6|) of the telegraph approximation is negligible. In this 
case the solution of the telegraph equation behaves like 
the solution of the diffusion equation. At long times the 
telegraph and diffusion approximations are almost undis- 
tinguishable and both are good, while the diffusion one 
is slightly better. 



FIG. 8. (Color online) The densities Punwci^^^) and 
PoRWBi^h^) and the telegraph and diffusion approximations 
PTE(a;,f) and p-Q-g^{x,t) at t — 20, a = 1, a — 0.1 and 9 — 1 
(e = 0.5). The telegraph and diffusion approximations are 
very close to each other. The singular terms Pdrwg and p^E 
are depicted by vertical arrows. Note that the singular term 
Pdrwg is quite "heavy". 
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FIG. 9. (Color online) The densities Pdrwg(^>^) and 
PnawBi^Wt) and the telegraph and diffusion approximations 
P'Y^{x,t) and pj3E(a;,t) a,t t — 50, t = 100 and t = 500, a — 1, 
a = 0.1 and e = 1 (e = 0.2, e = 0.1 and e = 0.02). The 
telegraph and diffusion approximations are almost undistin- 
guishable in the figure. The singular terms Pdrwg and pxE 
(only for t = 50) are depicted by vertical arrows. Cf. with 

Fig.ni 

VI. CONCLUSIVE REMARKS 

We should remind that the telegraph equation de- 
scribes transport phenomena with finite propagation 
speed, while the velocity of the motion of the particle in 
the DCTRWs is infinite, i.e., the DCTRWs do not simu- 
late the most distinctive property of the telegraph equa- 
tion. Moreover, asymptotic analysis and computations 
performed in this paper show that the diffusion equation 
gives better approximations to the DCTRWs than the 
telegraph equation. In other words, the DCTRWs are al- 
ways closer to their continuous limit than to the solutions 
of the telegraph equation. This implies that, in contrast 
to the widespread opinion, the DCTRWs do not simulate 
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the telegraph equation at the microscopic level. Likewise, 
CTRWs with the exponentially distributed and general 
waiting time do not simulate the telegraph equation at 
the microscopic level. 

Summing up, we can conclude that, first, the DCTRWs 
do not provide any correct microscopic interpretation of 
the one-dimensional telegraph equation, and second, the 
kinetic (exact) model of the telegraph equation is differ- 
ent from the model based on the DCTRWs. 

An interesting question arises concerning discrete- 
space random walks with the distributions of the 
jump length Ad, Eg. (IIV.4p . and the waiting time tptu 
Eq. (jIV.5[) . For sufficiently large /i this walk is very close 
to the Bernoulli random walk, Eq. (jIV.4|) . For /i = 1 
this is the discrete-space CTRW with the exponentially 
distributed waiting time (1/T)e*/'^. In the former case 
the telegraph approximation is better than the diffusion 
one, while in the latter case the diffusion approximation is 
better. Being incompletely stated yet, the problem is the 
following one: is there such that for fi e [l,A*o) the 
diffusion equation gives better approximation to these 
discrete-space CTRWs than the telegraph one, while for 
fi £ (/iOjOo) the telegraph approximation is better? 
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Appendix A: Two particular (discrete- and 
continuous-space) DCTRWs and the telegraph and 
diffusion approximations 

The density of the DCTRW with the distribution Ab 
(|V.1I) of the jump length results from the discrete-space 
variant of the formula pi.Sp . We denote the density by 
Pdrwb- It is given by 

CO 

PURWB{xi,t) = ^ Pa,nit)XB,nixi), 
n=\l\ 

xi =la, l£ Z, (A.l) 

where the probability density Ae^n is given by the scaled 
binomial distribution 



= < 



2"[(n-t-|;|)/2]![(n-|/|)/2]!' 

|Z| < n and n -I- ^ is even, 

0, n < |Z| or n -I- / is odd. 



and pa,n stands for the probability p„, corresponding to 
the distribution ■ The Laplace transform of Pa,n is 



CpaAs) = (^) 



-2/9 



[(s + 1/61)2 -(l-a)/6'2]«+i 
Therefore, the probability pa,n is given by 



e-*/« f ^ 



fit) + -J{t) 



< a < 1, 



1 ft 



2n 



(2n)! \e 



(2n+l)! \e 



2n+l 



a = 1, 



where 



|^(— 1)™ gVi-f t/s _j_^_]^y«+l g-\/l-a t/t 



with corresponding limit as a — ^ 1. 

The density of the DCTRW with the Gaussian dis- 
tribution Aq (jV.2p of the jump length results from the 
formula pLSp . We denote the density by Pdyiwg- 
given by 

Pdrwg(^'0 — Pdrwg(^'0 + Pdrwg(2^' (^-2) 
where the singular term is 

PLRWG(a;,i) =f'„(t)5(x), (A.3) 

with 



-t/e 



t 



cosh ( — 



1 



and the regular term is 



sinh ( --\/l — a 



t 



< a < 1, 
a = 1, 

(A.4) 



pWwoM = 5:Pa,„(t)^=^e-V(2-^) 

V27rncr 

n— 1 ^ 

The telegraph approximation is given by the solution to 
the initial- value problem (jIIL4p . (jIIL5|) for the telegraph 
equation 0, [H, SI] 

PTE,ix, t) = Pt^{x, t) + p^^{x, t), (A.5) 

where the singular term is 



p^rsix, t) = e-*/^ - [5{x - vt) + S{x + vt)] , (A.6) 
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and the regular term is 




(A.7) 



V — \faa I is the velocity, H{-) is the Heaviside 
step function. 

The diffusion approximation is given by the solution to 
the initial- value problem (|IV.2p . (|IV.3p for the diffusion 
equation 



(A., 



Note that the significant difference between the density 
Pdrwg of the DCTRW and the solution p^E (EH) 

of the telegraph equation is that the support of Pdr^vg 
(jA.3P is localized at the starting point of the walk {x — 
0), while the support of p^g (|A.6p is localized at the 
moving front (|a;| = vt). The term PDf{;wG ^'^ negligible 
for at/9 1, while the term is negligible for t/^ ^ 1, 
the former condition being stronger. Note also that the 
support of /Odrwg('' ^) '^hile the support of Pxe('' 
is [—vt,vt]. However, the latter difference appears to be 
in general not so significant. Indeed, according to the 
DeMoivre-Laplace and the central limit theorems psl-ls^] 
the solution (IA.8|) of the diffusion equation approximates 
the Bernoulli random walk (jIV.4l) with h = ^/aa and 
r = 29, t/{29) in the expression (|A.8[) being the number 
of jumps. The solution (jA.8| has the support R, while 
the particle involved into the random walk has a finite 
velocity, nevertheless, the diffusion approximation is a 
classic and widely used one to the Bernoulli random walk. 



at long times (the integrals are taken over the interval 
(0,00), rather than (—00,00), since the integrands are 
even with respect to k) . 

The Fourier-Laplace transform of the density P^j^wq 
of the DCTRW with the Gaussian distribution Aq of the 
jump length is the following one 

•^■^PdrwgI^' 

s + 2/9 



{s + 1/61)2 _ [1 _ „(i _ e-'^=fcV2)]/6i2 ' 



which results directly from the Montroll-Weiss formula 
pi.2l) with J-Xq and C^a ■ This implies the Fourier trans- 
form 



DRWG 



ik,t) 



-t/e 



cosh — ar(fc)^ 



1 



r{k) = 1 - e-'^''='/2, 
the Fourier transform of the singular term being 

-^PDRWG(fc:i)=^'a(i), (B.2) 



where is given by Eq. (IA.4[) . Obviously 

•^Pi3rwg(^' = •^Pdrwg(^'*) ~ •^PdrwgC''' ^) — >■ as 
/c — >■ 00. 

The Fourier-Laplace transform (jin.6|) of the telegraph 
approximation can be rewritten as 



2/61 



(s-f 1/61)2 -(l-a(72fc2/2)/6'2 ■ 
This implies the Fourier transform 



Appendix B: Asymptotic comparison of the 
telegraph and diffusion approximations to the 
continuous-space DCTRW at long times 

In this section the long time means at/9 ^ 1, i. e., time 
is longer than the mean waiting time 29/a. We com- 
pare the telegraph and diffusion approximations to the 
continuous-space DCTRW at long times by evaluation 
the L2(R) norms of the differences Pte('' ~ Pdrwg 
and Pde(''^) ~ Pi3RWg(''0- It is more convenient to 
consider the L2(K) norms of the Fourier transforms of 
the differences. Recall that the -L2(K) norms of a func- 
tion / and its Fourier transform J- f are related by 
= 11-7^/11 ^^(R)- Thus, we need to estimate 
the integrals 
1*00 

I I-^Pte(^:0 - •^PDRWG(*'Ol^dfc and 
Jo 

POO 

/ l-^PDE(fc,i)--^PDRWG(^-Ol'dfc (B.l) 

Jo 



-t/e 



cosh I — a—^— 




1 -a- 



sinh I -\ 1 



cr^k^ 



where we assume that if aa k > 2, then 
■\/l — aa-^k^/2 — iy^ aa^k^ /2 — 1. The Fourier transform 
of the singular term is 



Tpf^Eik^t) = e-*/^cos(wtfc) 



-t/e, 



t 'Jaak 
9 V2 



(B.3) 



Obviously Tp\;^{k,t) = J^pj.^{k,t) — J'p?p^{k,t) — as 
k — 00. 

The Fourier transform of the diffusion approximation 

is 

j-pj^j^(fc,t) = e-('^*/^)(-''='/4). 
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The "singular terms" J^Py^^^j^/q (|B.2[) and J^p^E 
are exponentially small as t — >■ oo and, therefore, neg- 
ligible on any finite interval with respect to k for large 
t. 

To estimate the integrals (jB.ip at long times it is con- 
venient to introduce the small parameter e = (at/ 9)^^ <C 
1 and change the variable k for C, = ak/{2^) — 
{ak/2)y/^aje. Then we have in the new variable the 
Fourier transform of the density of the continuous-space 
DCTRW 



'drwg(C' i) — e "^/"^ 



cosh 



ae 



: sinh 



v/1 - ar{C) 



ae 



r(C) = 1 



-2eC 



being ■FPdr wg(C.O 
sinh (Vl — a/ ae) / \/T 



the "singular term" 

g-i/ae |^(,Qgj^ (VI - a/ae) + smn [V ^ - a/ae)/\/i - a \ 
the Fourier transform of the telegraph approximation 



as C — ^ oo- The reason for such a behaviour is that the 
regular term p^^{x,t) (|A.7p is discontinuous at |a;| — vt. 
Nevertheless, due to exponential decrease of ^Pte(Ci^) 
for small e the first integral of the integrals (|B.5p can 
be approximated with arbitrary accuracy by the integral 
with the upper limit instead of oo, if is sufficiently 
large. Note that for any Co there exists sufhciently small 
e such that Co < l/-\/2a£, i.e., [0, Co] C [0, l/\/2ae\. 

Thus, for small e the integrals (jB.sp can be approxi- 
mated with arbitrary accuracy by the integrals with the 
upper limit Co instead of oo, where Co is sufficiently large. 

The straightforward calculations imply uniform 

asymptotics on the interval [0,Co]: 

= [l+[l{e-Q')+C'\e}o-^' +0{e^) 

as e — > 



and 



•^PTE(C,i) = e 



cosh 



Vl-SaeC^' 



ae 



yjl - 2aeC 



: sinh 



Vl - 2aeC^ 



ae 



(B.4) 



the "singular term" being -^PteCCiO = 
g-i/ae (^y/2/ae C), and the Fourier transform of 
the diffusion approximation 

J'PBEiCt) = c"'-' • 

The integrals (jB.ip . within a factor of 2/(t, become in the 
new variable 

POO 

Ve l-^/3TE(C,i) --^PDRWolCOI^dC and 
Jo 

poo 

Vi / l-^PDE(C,i)--^PlDRWG(C,i)l'dC. (B.5) 

^0 

The "regular term" J"Pdrvvg(C,0 = -T^PdrwoIC,*) - 
•^/'drwg(C' ^) tends to zero monotonously and expo- 
nentially as C — ^ oo, as well as the Fourier transform 
J^Pj3g(C, i). Therefore, the second integral of the inte- 
grals (IB. 51) can be approximated with arbitrary accuracy 
by the integral with the upper limit Co instead of oo, if 
Co is sufficiently large. 

For small e the Fourier transform -^Pxe(Ci^) tends 
to zero monotonously and exponentially on the interval 
[0, l/^/2ae], where the square root in Eq. (|B.4[) is nonneg- 
ative. The "singular term" J^PxE '^^'^ be neglected on this 
interval for small £. Asymptotic behaviour of the "regular 
term" is J^Pte(C, t) = -^PteIC, - •^'PteIC, t) = 0(C"^) 



^p!pE(C,i) = {l+g {e-C')\e]^~'^+0{e') 



as e 



(-^/'drwg(Cj J^prpg(C, have the same asymp- 
totics). Therefore, the asymptotics for the differences 
are 



•^PTE(C,i) - -T^PdrwgIC,*) 

= -C^e-'^'e-|-0(£2) as e^O (B.i 

and 

•^Pde(C, i) - -7^Pdrwg(C, t) 



^ (C' - C*) + C'J e-^ e + 0(e2) as e ^ 0. 

(B.7) 



Therefore, at long times (for small e) the integrals 
(lR5t are 



POO 

V~e |J-p^E(C,i)--^p|DRWG(C,i)l'dC 
Jo 



.2-5 / c«c-2C^dC = yi^e^-^ 







Q.257e^'' (B 
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and 



10 



-2.5 



dC 



3!! 5!! 7!! 

27 ~ 28 2" 



5!! 7!! 



7!! 



-2.5 



2^ 2V 29 
(0.020a2 - 0.110a + 0.257) e^ ^ (B.9) 



Note that at long times the same integral for the differ- 
ence ^Pte ~ -^P 



IS 



/•oo 

l-^PTE(C,i)--^PDE(C,Ol'dC 

Jo 



-2.5 



-2C^ 



dC « 0.020 £2 



(B.IO) 



The asymptotic estimates (jB.Sp and (|B.9p show that 
at long times the diffusion approximation to the density 
Pbrwg '^^ DCTRW is better in the 1^2 (K) norm, than 
the telegraph approximation. At the same time, for small 
a the telegraph approximation is almost as good as the 
diffusion one. 
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